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Abstract 



Real-time perturbation theory is formulated for complex scalar fields away 
from thermal equilibrium in such a way that dissipative effects arising from 
the absorptive parts of loop diagrams are approximately resummed into the 
unperturbed propagators. Low order calculations of physical quantities then 
involve quasiparticle occupation numbers which evolve with the changing state 
of the field system, in contrast to standard perturbation theory, where these 
occupation numbers are frozen at their initial values. The evolution equa- 
tion of the occupation numbers can be cast approximately in the form of a 
Boltzmann equation. Particular attention is given to the effects of a non-zero 
chemical potential, and it is found that the thermal masses and decay widths 
of quasiparticle modes are different for particles and antiparticles. 
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I. INTRODUCTION 



Many problems in physics require an understanding of the non-equilibrium properties 
of highly excited states of a quantum field theory. The principal motivation underlying the 
present work is to understand the nature of phase transitions in the very early universe, in 
particular those that might give rise to inflation |I|-§]|. Most discussions of inflation assume 
that the non-equilibrium state of the relevant system of quantum fields can adequately be 
characterised by expectation values governed by a classical potential. However, calculations 
which attempt to go beyond this classical approximation indicate that it may be seriously 
inaccurate (see e.g. Non-equilibrium effects are also likely to be important in the 

context of baryon-number violation at the electroweak scale |9|,|l(J and in the study of quark- 
gluon plasmas formed in heavy- ion collisions [II|, as well as in many aspects of condensed 



matter physics. 

Our own long-term programme, set out in \V2], has as its goal a complete description 



within perturbation theory of the dynamics of symmetry-breaking phase transitions. An 
essential requirement is to formulate perturbation theory in such a way that the evolving 
state of the non-equilibrium quantum fields is taken adequately into account in the low orders 
of the perturbative expansion which are likely to be tractable. The standard closed-time- 
path formalism flT5|-fnj|1 which deals with real-time phenomena in thermal equilibrium was 
generalised by Semenoff and Weiss [19 to apply to a real scalar field with a time- dependent 
action, such as might arise in a Robert son- Walker spacetime. However, the propagators 
entering their Feynman rules depend on the state of the field essentially through single- 
particle occupation numbers which are fixed at their initial values, and therefore do not 
properly reflect the evolution of the state. 

In earlier work [20)- 22] we have shown, for real scalar fields how this difficulty may be 
overcome. By adding a suitable counterterm (which we call the "dissipative" counterterm) to 
the quadratic part of the action, and subtracting it from the interaction, we obtain a lowest- 
order theory in which the absorptive parts of loop contributions to the 2-point functions are 
partly resummed into the unperturbed propagators. As a result, these propagators contain 
quasiparticle occupation numbers which evolve with time in the expected way. Interestingly, 
the form of the counterterm is rather similar to the action encountered by Hu and various 
collaborators (see e.g. PBH2H] and references given in these papers). These authors use the 



Feynman- Vernon influence functional technique to obtain the effective action for a quantum 
field in the presence of environmental variables (which can be thought of as a heat bath) 
when the latter are integrated out. Our procedure can be thought of as a self-consistent 
treatment of a field which provides its own heat bath, though none of the modes of this 
field are integrated out. A similar structure also arises in non-equilibrium formulations of 
thermo- field dynamics (see e.g. ||26|| ) but this formalism turns out to equivalent to a version 
of quantum statistical mechanics restricted to Gaussian initial density matrices |27| 

In dealing with a realistic model (such as the standard model or a grand unified theory), 
one naturally meets complex scalar, spinor and gauge fields and it is essential to extend 
the non-equilibrium formalism to encompass fields of these kinds. (Real-time perturbation 
theory for scalar, spinor and gauge fields in equilibrium is discussed, for example, by Kobes 
et al [f^J.) The purpose of the present work is to extend the formalism to complex scalar 
fields, giving particular attention to the effects of a non-zero chemical potential, which 
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turns out to be non-trivial. In section II we obtain the structure of the matrix of 2-point 
functions in the closed-time-path formalism, which is applied in section III to find the most 
general permissible form of the dissipative counterterm. The corresponding unperturbed 
propagators are derived in section IV, where we find that, in the presence of a non-zero 
chemical potential, not only occupation numbers but also decay widths and thermal masses 
are different for the quasiparticle modes corresponding to particles and antiparticles. In 
section V we evaluate the dissipative counterterm at the lowest non-trivial order and find, 
with suitable approximations, that the time-dependent occupation numbers obey a kinetic 
equation of the Boltzmann type. Finally, our results are summarised in section VI. 



II. PROPAGATORS IN THE CLOSED-TIME-PATH FORMALISM 

We consider the theory of a massive, complex, self- interacting scalar field, 0(x, t), de- 
scribed by the Lagrangian density 

£(0, 0*) = 0*0 - V0*V0 - m 2 (t)0*0 - ^(0*0) 2 (1) 

where = d<p/dt. The generating functional for the time ordered Green's functions of the 
theory may be written as 

Z\j] = Tr \ p Te i f d4x{j ' t '* +j *' t ' ) ] (2) 

where 0(x, t) is the Heisenberg-picture field operator, j(x.,t) is a complex source and T 
denotes time ordering. For an initial state of thermal equilibrium at time to, characterised 
by the temperature 1/(3 and chemical potential /i, the density operator is 

e -/3[H(t )-tiN] 

P = Tr [g-/3[ff(t )- M iV]] ' ( 3 ) 

where H(t) is the time-dependent Hamiltonian and N is the particle number operator. It 
is convenient to generalise this generating functional, using the Schrodinger picture field 
operator s (x) to write 

Z[ji,j 2 ] = Tr pTe i f d M#s+j20*+j*<My e -i J d 4 x(H s -ji^*-j^<f, s ) 

where T denotes anti-time ordering and H s (t) is the Hamiltonian expressed in terms of S . 
This Hamiltonian depends explicitly on time through the time-dependent mass m(t). The 
original generating functional (Q) is given by Z[j] = Z[j,0}. 

Each of the three factors inside the trace in (§) has a path-integral representation with 
its own integration variable, say 0i(x, t) for the time-ordered factor, 02 (x, t) for the anti- 
time-ordered factor and 03 (x, r) for the density operator, where the real time t runs from 
the initial time to to some large final time tf, while the imaginary time r runs from to 



(3. As usual, [p]9| , |20| , these three fields can be envisaged as living on the three segments 
of a contour in the complex time plane illustrated in [TJ. The effect of a non-zero chemical 
potential can formally be taken into account by using H{to) as the generator of evolution in 
imaginary time and imposing the boundary conditions 
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3 (x,/?) = e-^<Mx,t ), ^(x,/?) = e^#(x,t ) • 



(5) 



On introducing a third source for 3 , the generating functional @ is then given by the path 
integral 



Z\jij2j3]=^J[m[dM[^]e iS ^'* a) 
where M is a normalizing constant and 

2 , 3 ) = jf ' dt J d 3 x [£(00 - £(0 2 ) + j 1( f>l + fifa + j 2 <j>% + fM 



(6) 



+ J Q drJ d 3 x [i£ E (<fa)+h<i>$+jM 



(?) 



Here, Ce denotes the euclidean version of (|J). 

In order to carry out perturbation theory, we wish to split the action into an unperturbed 
part, Sq, an interaction part Sj and a source term Sji 



S((f)l, 02, 03) = 5'o(01, 02, 03) + 5/(01, 02, 03) + Sj(0i, 2 , 3 ) • 

In the standard way, the generating functional (§) can then be rewritten as 



Zo[ji,j2,j3] = exp 



d xd y j*(x)g ab (x } y)j b (y) 



The unperturbed propagator g(x,y) satisfies 

V k (x, d x )g{x, y) = -i6(x -y) = g(x, y)V(y, -d~ y ) 



(8) 



(9) 
(10) 



(11) 



where T>(x, d) is the differential operator corresponding to Sq, i.e. So = 

/ c / 4 x0:(x)p a6 (x,9)0 6 (x). 

The perturbation theory which results from choosing Sq to be simply the quadratic part 
of S was derived some time ago by Semenoff and Weiss [19| for a theory with somewhat 
more general time dependence than (|l|). It has the disadvantage that coefficients in g(x,y) 
which correspond roughly to particle occupation numbers are fixed at their initial values, 
and do not change so as to reflect the evolving state of the system. As shown by one of us 
for the case of a real scalar field 0,^1 > this situation can be improved by choosing instead 



S O (01,02,03) = S (2) (01, 02, 3 ) + / d 4 X<P* a M ab {x,d)^ 



(12) 



where M is a differential operator, to be chosen in such a way that the unperturbed prop- 
agator g(x,y) mimics as nearly as possible the dissipative behaviour of the full propagator 
G(x, y). To make perturbation theory tractable, we require the term involving M to be local, 
so it can involve only the real-time fields 0i and 02- Since (as is worth emphasising) we do 
not wish to change the overall theory, the interaction term Si now includes a counterterm 
—<p* a M ab (j) h in addition to the original interaction — (A/4)[(0];0i) 2 — (0202) 2 ] and M will be 
chosen so that this counterterm subtracts some part of the loop contributions to the full 
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propagator. In this way, we partially resum the absorptive parts of these loop contributions, 
and optimise g{x,y) as an approximation to G(x,y). In particular, we will find that g{x,y) 
now involves quasiparticle occupation numbers which evolve with time in the expected way. 

In order to determine the permissible form of M, we first investigate the structure of the 
full propagator, whose real-time components (a, b = 1,2) are given by 



Gab X J 











j=0 



dj*(x) dj b (x' 
(T[<f>(x,t)<f>*(x',t')})„ (0*(x',O0(M)>m 

(0(x,t)0*(x',t')>M (f[0(x,t)0*(x',t')]>M 



(13) 



where (• • -) M indicates a thermal average in the ensemble with chemical potential /i. Writing 

(0(x)0*(x')) M = n„(x, x')Q(t - t') + K„(x, x')Q(t' - t) (14) 

it is simple to show from the Hermiticity of p that /C M (x, x') = TC*(x', x) and from the charge 
conjugation (C) symmetry of ([!]) that 



(4>*(x')<f){x)) ll = H*(x, x')Q(t - t') + H-^x', x)G(t' - t). 



(15) 



For any complex function of fi, we define the Jj-conjugate /"(A*) = /*( — A 4 ) 5 which is related 
to charge conjugation. With this notation, the full propagator can be expressed as 



G(x,x') 



G(t - if) 



H(x,x') n ( (x,x') 
H(x,x') H\x : x') 

H*\x',x) Ti*\x',x) 
H*(x',x) H*{x',x) 



+ 



Q(t' - 1) 



(16) 



where the suffix /i has been suppressed. This structure should be maintained if the scalar 
field is embedded in a more general C- or CP-invariant theory. 



III. THE DISSIPATIVE COUNTERTERM 

We now want to construct a differential operator T> and the corresponding counterterm 
matrix M in such a way that equation (|TT|) admits a solution for g(x,y) of the form (|1~6|). 
Since we are dealing with a spatially homogeneous theory, we can take a spatial Fourier 
transform, obtaining a propagator gk(t,t'), and operator Vk(t,d t ) which obey (|TT|) in the 
form 

T> k (t, dl)g k (t, t') = -i5(t - = g k (t, t')V k {t\ -%) ■ (17) 

In the following, we will usually suppress the suffix k. As in |20|| , we first construct the 
operator T>(d t ) with constant coefficients appropriate to a temporally homogeneous system, 
and then allow these coefficients to depend on time. In the temporally homogeneous case, 
g(t,t') can be expressed in the form of (p"6"|) in terms of a function h(t — t') (which also 
depends on k and /z) for which we write 
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h(t-t') = u(t-t')+iv(t-t') (18) 
h*$(t-t') = w(t-t')+iz(t-t') (19) 

where u, v, w, and z are real functions. (Note that /*"(//) = /(—A 4 ) for any function /(//) ). 
For the temporal Fourier transform 

/oo 
dte^git) (20) 
-oo 

we obtain 

gn{u) = g 22 (u) 

= [A{u) - b(u) + C{u) + d{u)] 

+i [a(u) + B(u) - c(u) + D(u)] (21) 
g 12 (u) = 2[C(u) + d(u)] 
g 21 (u) = 2[A(u)-b(u)] 



where 



roc 

A(oj) + iaiuj) = \ dtu(t) [cos tut + ism Lut] 
Jo 

POO 

B(uj) + ib{uo) = / dtv(t) [cos ut + i sin ut] 
Jo 

roo 

C(u)+ic(u)= / dtw(t) [cos ut + i smut] 
Jo 

D(u) + id(u) = / dtz(t) [cos cut + i smut] 
Jo 



Note that the change of variables {u,fi) — > (—u,—fi) leads to the interchange 
(A(u), a(u), B(u), b(u)) <-> (C(u), —c(u), D(u), —d(u)), and consequently 



g ab (u) = g* b l(-u) (22) 



which also follows directly from (0). Upon Fourier transformation, ([17|) becomes an alge- 
braic equation whose solution for T>(—iu) is 



V(-iu) 



det \g\ 



922{U) -g\2{0J) 
-hi (w) (a;) 



(23) 



This solution shows how T> can be constructed from the (unknown) functions A(u) ■ ■ ■ D(u) 
which are real, even functions of u and a(u) ■ ■ ■ d(u) which are real and odd. We will 
choose these eight functions in such a way that T> has certain essential properties. First, 
we require T>(dt) to be a second-order differential operator. The 5(t — t') in ([17]) arises 
from the derivatives of Q(t — t') and Q(t' — t) and, in order that these appear only in the 
diagonal elements, second derivatives may appear only in the diagonal elements of T>(d t ). 
To maintain the normalization of (f) as in (P, we fix the coefficients of <9 t 2 , or of — u 2 , to be 
±1. Expanding ( [23] ) in powers of u, we find that the most general form of T> satisfying these 
requirements can be expressed in terms of six real coefficients a, /3, 7, a, 7, 7 as 
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-uj 1 — 7(7 — i^y)u + (3 



i a 



-z(7 — 7)cj + z(a + a) 



z(7 + 7)cj + i(a — a) c<j 2 — 2(7 + i^y)u — [3 — ia 



(24) 



Since V(—iu) is equal to — ig 1 (u>), it has the same symmetry fl22|) as g(u;), and we may 



deduce that a, /3 and 7 are even functions of /z, while a, 7 and 7 are odd. Since these 
quantities are also real, we have 



a" = a, /3 B = 0, 7 B = 7 
a" = —a, 7" = —7, 7^ = —7 



(25) 



In particular, when fi vanishes, we have a = 7 = 7 = and P reduces to the form found in 
|2(J for a real scalar field. 

For the non-equilibrium theory, with time- dependent bare mass, T>(t, dt) is given by 



T>u = d 2 t + (7 - ij)d t + p - ia + (7 - i7)/2 
£>i 2 = (7 - 7)^ + i(Q! + a) + (7 - 7)/2 
^21 = -(7 + 7)$ + i(a - a) - (7 + 7)/2 
^22 = + (7 + 17)$ - /3 - + (7 + fy)/2 



(26) 



where the coefficients are now time dependent, and the derivatives 7, etc ensure that T>(t, dt) 
is a symmetrical operator, as required in (JTTJ) . That is, jd t + 7/2 = 7 1 / 2 <9t7 1 / 2 , etc. We 
emphasise again that the six coefficients in (p6|) , while undetermined at this stage, do not 
represent arbitrary modifications of the theory, but rather an optimal choice of the lowest- 
order theory which is to serve as the basis for perturbation theory. As indicated above, these 
coefficients will be found self-consistently from a suitable renormalization prescription, and 
we shall find that they correspond to six properties of a gas of quasiparticles, namely the 
thermal masses, occupation numbers and relaxation rates of particles and antiparticles. 



IV. THE DISSIPATIVE PROPAGATOR 



The basic equation (|i~7D for the unperturbed propagator matrix g(t,t') is satisfied if the 
function h(t, t') obeys the conditions 



d t [h(t,t')-h**(t',t)] t ,_ t 

dt[h(t,t')-h*(t',t)]\ t , =t =0 

[V u (t,dt)+V 12 (t,d t )]h(t,t') = 

V n (t, d t )h*\t', t) + V^t, d t )h\t', t) = 

together with their jj-conjugates and the continuity conditions 

h(t,t) = h*{t,t) = h\t,t) . 

We find that the solution of these equations can be written in the form 



(27) 
(28) 
(29) 
(30) 

(31) 



h(t,t' 



4:Jtt(t)tt(t') 



(l + W(0)e -i -# n(i) * 



(32) 
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in terms of the auxiliary functions Q(t) and N(t) which are solutions of 



in 3n 2 o2 o i. . _. 2 . 



(33) 



d t + 2%n - ~ + 7 



with the subsidiary conditions 



[(<9 t + 7)iV - 7] = 2za + 2iaiV + ry(7iV - 7) 



(34) 



n 



n 



(35) 




1 /n n* \ fN + N^ 

-2 U + sf +27 hr- 



— z 



(iV-iV tt ) - (N-N*)' 



(36) 



It is straightforward, though tedious, to show that these conditions are preserved by (|34]) if 
they hold at the initial time. 

With a non-zero chemical potential, 7 and a are in general nonzero, so the frequency 
n(t) which satisfies (|33]) is complex, though it has the property ^ = Q. 

), we require the full 3x3 



To obtain the initial conditions that apply to (33) and 
matrix of propagators which satisfies an equation of the form QT7D with 



(37) 



and £>i 3 = P 2 3 = ^31 = ^32 = 0, where u 2 = k 2 + m 2 (t ). The boundary conditions 
which apply to these nine propagators are set out in the Appendix, where results for those 
involving imaginary times are also given. 

For the case of a time- independent mass m 2 , we expect a steady state solution, with 
Q = N = 0. In this case, we find that the boundary conditions can be satisfied only if three 
relations hold between a, (3, • • • . The first of these is 



f 



A{(3-u 2 ) . 



(38) 



It seems natural to require that the renormalized masses (or, more generally, the fc-dependent 
frequencies) in T>\\ and D33 should be equal, although this is not obligatory. This means 
that P = uj 2 = m 2 + k 2 , and hence that 7 = 0. We will indeed assume that 7 = 0, since this 
introduces considerable simplifications, but the consequences of this assumption will need 
some discussion at a later stage. In this case, the two other relations are 



a + a — (7 — 7)0;] — (70? — at)n 



1 
2 

- [a — a — (7 + 7)0;] — (70; + a)n H 







(39) 
(40) 



where 
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are the usual occupation numbers for particles and antiparticles in a free Bose-Einstein gas. 
The coefficient N in the real-time propagators is then given by 



N= f ° + ' 7/2 I (l+n+ + n -) + (n- - n + ) (42) 



and the frequency fl by 



Q 2 = oj 2 - 774 + ia . (43) 



Note that when the dissipative counterterm is neglected, we have Q = uj and 7 = 0, so that 
N = 1 + 2rT . We then find, as expected, that the coefficient of the positive frequency term 
in the propagator (^) is (1 + N) = 2(1 + n~) while that of the negative frequency term 



is (—1 + N^) = 2n + . Thus, with the dissipative terms present, the real-time propagators 
describe a gas of quasiparticles, in which quasiparticle modes contain a small admixture of 
bare antiparticles and quasi-antiparticle modes contain a small admixture of bare particles. 

For this gas of quasiparticles, the positive- and negative-frequency mode functions are 
exp[— iQ + (t — t')] and exp[iQ_(t — t')] respectively, where 

fi± = [Reft =p 7/2] =F i [7/2 t Imfi] . (44) 

We see that the decay rates for particle and antiparticle modes will be different if Imfl ^ 0, 
which will in general be true if either 7 or a is nonzero. On the other hand, the thermal 
masses of quasiparticles and quasi-antiparticles will be different if and only if 7 is nonzero. 
Which of these conditions actually applies will be discussed in the next section. 

For the non-equilibrium theory, any solution of fl33|) may in principle be used for Q(t). 
Clearly, however, it is desirable that our lowest-order theory should approximately retain 
the characteristics of the steady-state solution in the case where m 2 (t) is slowly varying near 
to. We will therefore choose a renormalization prescription for which /?(t ) — and the 
solution of ( p3|) which has Cl(t Q ) = and Q (t ) given by (f|3|). Retaining the assumption that 
j(t ) = 0, (^) shows that Cl(t ) = also. We would now like to interpret N(t) in terms of 
time-dependent occupation numbers n ± (t). To this end, we define 

nHt) = {wm ± n(t)mmt) + 

Ttl{t)[N{t) - N\t)] - 2Q{t)} (45) 

so that 

N(t) = ( m ^ )/2 ) [l + n + {t) + n-{t)] + [n-{t) - n + (t)} . (46) 

It is easy to see that n + (t) and n~(t) are real and that n + ^(t) = n~(t) as required for the 
occupation numbers. If n ± (t ) are taken to be the equilibrium values (f|l]), then the desired 
initial condition N(t ) = ensures that (|36|) is satisfied. 
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As in [20], we can now show that with reasonable approximations, the evolution equation 
(|34l) and its (^-conjugate reduce to a pair of Boltzmann-like kinetic equations for n (t). We 
first rewrite (|3~4|) with 7 = as 



Q [d t + 2i(Q - 17/2)] P = 2i(a + a)N 

p = n- 1 [(d t + 1 )N-^} 



(47) 
(48) 



and assume that dtP « (2Q — ij)P, which will be valid if the characteristic relaxation 
time I/7 is much greater than 1/Q. It will turn out, as in [p0~| |, that 7 is smaller than Q by 
a factor of A 2 . Indeed, 7, 7, a and a are all of order A 2 . Consequently, to order A 2 , we can 
replace N by its lowest-order value N rs (1 + 2n~) and neglect 27/2 in comparison with Q. 
With these approximations, ([|7]) and its jj-conjugate become 



dn 


1 


dt 


2ft 


dn + 


1 r 


dt 


= 2n [( 



[a + a - (7 - 7)^] - ^(7^ - a ) n ~ 
a ~ (7 + 7)^] - ^(7^ + «)™ + 



(49) 
(50) 



Evidently, the relations (|39|) and (]40[) which apply to the steady state solution are just the 
conditions for the occupation numbers to be constant. In the next section, we evaluate a, 
a, 7 and 7 explicitly, and find that the right hand sides of (|49|) and (|50D have approximately 
the form of the scattering integrals which appear in the Boltzmann equation. 



V. EVALUATION OF THE DISSIPATIVE COUNTERTERM 

As explained above, our strategy is to choose the counterterm M so as to optimise the 
unperturbed propagator g(x,y) as an approximation to the full propagator G(x,y). The 
relation between these may be expressed by the Dyson-Schwinger equation 

G ab (x,y) = g ab (x,y) +i d 4 zd^z' g ac (x, z)i: cd (z } z')G db (z' ,y) (51) 



in terms of the self-energy y). As explained in detail in [20[, we express S in terms of 
the average t and difference At of its time arguments, and take the Fourier transform with 
respect to (x — y) and At. The result can be written as 

= M(t, -iu) + E«(f, u 2 )u + t^{i,u 2 ) (52) 

where S denotes the contribution from loop diagrams, which has been split into parts even 
and odd in the transform variable u. In the counterterm contribution, the i dependence is 
that of a{t) ■ ■ ■ while — iu replaces dt. While M(i, —iu) is a second-order polynomial in u, 
the same is not true of S(t, u). As in standard renormalization theory, therefore, we choose 
the counterterm to cancel E at some reference frequency, which can conveniently be chosen 
as fl = Q(t). This gives 

{ 1 (t)-m = ^ 2 \t,n 2 ) (53) 

Ht) + a(t)]=*Ei?M 2 ) (54) 
m-zm=^(t,n 2 ) (55) 
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from which a, a, 7, 7 and 7 can be deduced by isolating parts which are even and odd under 
jj- conjugation. The corresponding equation for /3, namely 



k 2 + m 2 (t) - /?(*) = ReX { 3(t, Q 



(56) 



can be adjusted to meet the requirement (3{to) = k 2 + m 2 R , where rriR is an appropriate 
renormalized mass, as discussed in the last section, but the details are not important for our 
present purposes. 

The lowest-order contribution to E 12 is from the graph shown in 0. An approximate 



method for evaluating it is described in [20] . Since the whole diagram is proportional to A 2 , 
we take the lowest-order approximation to its internal propagators, using 



h(t,t') 



2Q(t) 



(1 



+ n (t 



+ n + (i)e 



(57) 



where n ± (t) and Q(t) are assumed to vary sufficiently slowly that they can be treated as 



effectively constant. After calculating ak(t), ak(t), 7fc(t) and 7(2) from ( |53| ) 
substituting the results in ( |49"|) and (|50|) , we obtain the Boltzmann equation 



and 



d 
dl 



32(2tt) 



J d 3 k 1 d 3 k 2 d 3 k 3 5(ki + k 2 - k 3 - k) 



x8(Qi + fi 2 - fi 3 - [fiifi 2 ^3^ 



1-1 



x 



ntrq (1 + n+)(l + rq) - (1 + n+)(l + n 2 )n+n^ 
+n7n+(l + n+)(l + nj^ ) - (1 + n7/)(l + n+)n+ 
+% nj(l + ng )(1 + % ) - (1 + 7i7/)(l + nj)n3 



(58) 



where fii denotes Q kl (t), etc. At the order of approximation we are using, f2 fc (t) is real. 
The scattering integral on the right correctly describes the rate of production minus the 
rate of absorption of particles of momentum k due to all 2-body processes allowed by charge 
conservation. It vanishes when the occupation numbers have the Bose-Einstein form 



1 



n 



(59) 



where the chemical potential /2 need not be the same as the original //. Of course, the rate 
of change of n + is given by the Jj-conjugate of (|58|). 

The value of 7 obtained from ([)5|) is not zero although, unlike a, a, 7 and 7 h involves 
off-shell processes. As discussed in the last section, this implies that the thermal masses of 
particles and antiparticles are different when the chemical potential is nonzero, the difference 
being comparable in magnitude to the decay widths (7/2±Imf2) of the quasiparticle modes. 
We are nevertheless entitled to set 7 = 0. This simply means that the difference in thermal 
masses is treated purely perturbatively, rather that being resummed into the unperturbed 
propagator. As it happens, the contributions involving 7 to both ([33]) and ( |34D are of order 
A 4 . Thus, in low order calculations, it would be consistent to retain a nonzero value of 
7 in (|32|), while setting it to zero in the subsidiary calculations needed to determine Q(t) 
and N(t). The approximations needed to obtain the Boltzmann equation in any case entail 
setting 7 = 0. 
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VI. SUMMARY 



We have extended to the case of complex scalar fields a formulation of non-equilibrium 
perturbation theory which partially resums the dissipative effects of loop diagrams into the 
unperturbed propagator. As a result of this resummation, low-order calculations using our 
modified Feynman rules reflect the evolving state of the non-equilibrium system through 
terms in the propagator which can roughly be interpreted in terms of quasiparticle occupa- 
tion numbers. As in earlier work, we find that the evolution equation for these occupation 
numbers can be cast approximately in the form of a Boltzmann equation. Kinetic equations 
of a similar form can, of course, be derived by other methods (see, for example |2D|r3T|) if 
one sets about finding the rate of change of occupation numbers defined in some appropriate 
manner. In an interacting theory, however, the definition of single-particle modes is gener- 
ally ambiguous. We, on the other hand, have set about calculating unambiguously defined 
Green's functions in an optimal manner. For us, the quasiparticle modes and their occu- 
pation numbers simply provide a convenient way of thinking about quantities which arise 
naturally in the course of these calculations, and the Boltzmann equation is a low-order 
approximation to a rather more complicated second-order evolution equation. 

The novel aspects of the work reported here concern the effects of a non-zero chemical 
potential. To the extent that the quasiparticle picture is valid, we find that not only the 
occupation numbers but also the decay widths and thermal masses of the quasiparticles 
are different for particle and antiparticle modes. In particular, the difference in thermal 
masses is a two-loop effect (in the A(0*0) 2 theory), which is not readily apparent in the 
usual perturbative treatment (see e.g. j31|), although it should be derivable from a two- loop 
calculation of the self-energy, even in equilibrium. 



ACKNOWLEDGMENTS 

D.B. would like to thank the University of Leeds for financial support in the form of the 
William Wright Smith research scholarship. 



APPENDIX A: 

Boundary conditions on the nine propagators g a b arise from two sources. First, boundary 
conditions on the fields (H) together with 0i(x, tf) = 2 ( x , £/), 02( x , ^o) — 03( X ; O) and the 
complex conjugates of these lead to 

9la(tf,t') = g 2a (tf,t') 9al(t,tf) = 9a2(t,tf) 

92a(t ,t') = 03a(O,tf) <?a 2 (Mo) = SasM) (Al) 

9ia(t ,t') = e^g 3a {0,t') <? a i(Mo) = e-^ a3 (t,0) 



for (a = 1,2,3). Second, the evaluation of the Gaussian path integral leading to (|T0|) is 
performed, as usual, by completing the square and this involves integrations by parts. The 
requirement that boundary terms arising from these integrations by parts cancel leads to 
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dt9ia(t,t')\ t=tf = d t g 2a (t,t') 

1 „ 1 

2 



t=t f 



dt92a(t,t')\ t=to + -{ 1 + j)g la (ta,t') - -(7 + ij)g 2a (to, t') 



id T g3a(T, £')lr=0 + ^ai®, t') 



1 1 

d t gia(t, t') \ t=to + - (7 - -f)g 2a (to, t) + - (7 - ij)gia(to, t ) 



id T g 3a (r, t')\ T=/3 + -g 3a ((3, t') 



(A2) 



(A3) 



(A4) 



for (a = 1, 2, 3), together with the Jj-conjugates of these. The coefficient 6 corresponds to an 
additional counterterm which, as explained in detail in 21], enables us to impose N(t ) = 
as befits an initial equilibrium state. As in fl2"I |, we find 5 = —7. Here, 7, 7 an d 7 are all 
evaluated at t . 

The real-time propagators satisfying these boundary conditions are described in the text. 
Of those involving imaginary times, three are found to be 



y) = i.{[n- + e(r-0 



and 



933 V > 



9i3(t,r) = g 23 (t,T) 

-\ r n (7(t)-i7©)d* 



n 



+ 6(r / -r)] e" (T - T ' } } 



(A5) 



4yJn(t )Q(t) 
x {A x e~^ 



A 2 ene~ lJt o nm +(A 3 e 



with 



A x = (1 + n + ) [iV(t ) - 1 - 2n~]/(l + n + + n~) 
A 2 = 7i~[iV(to) + 1 + 2n + ]/(l + n + + n~) 
A3 = (l + n + )[iV(to) + l + 2n-]/(l 
A 4 = n~[A^(to) - 1 - 2n + ]/(l + n + + rT) 



n + + n 



The others are given by 

93i(r,t) = g 32 (r,t) = g$(t,r) 
Indeed, the whole propagator matrix has the property 

9*l(t,t')=g ba (t,t) 

where t and t' denote real or imaginary times, depending on the indices a and b. 



(A6) 



(A7) 
(A8) 
(A9) 
(A10) 



(All) 



(A12) 
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FIGURES 

FIG. 1. The closed time path in the complex time plane 
FIG. 2. Lowest-order Feynman diagram contributing to the dissipative counterterm 
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